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^ , Abstract. We present the results of observations from the Plateau de Bure IRAM interferometric array of the 

submillimetre (submm) galaxy SMM J16359+6612 lying at z = 2.516 behind the core of the massive cluster 
A 2218. The foreground gravitational lens produces three images with a total magnification of 45 of this faint 
submm galaxy, which has an intrinsic submm flux of just /sso/jm = 0.8 mjy - placing it below the confusion 
' limit of blank-fleld surveys. The substantial magniflcation provides a rare opportunity to probe the nature of a 

I distant sub-mJy submm-selected galaxy, part of the population which produces the bulk of the cosmic far-infrared 

. background at submm wavelengths. Our observations detect the CO(3-2) line in all three images, as well as the 

CO(7-6) line and the dust continuum at 1.3mm for the brightest image but only at a 3o" level. The velocity proflle 
of the CO(3-2) line displays a double-peak profile which is well fit by two Gaussians with FWHM of 220kms~^ 
and separated by 280kms~^. We estimate the dynamical mass of the system to be ~ 1.5 x 10^" Mq and an H2 
O ' gas mass of 4.5 x 10^ M©. We identify a spatial offset of ~ 1" between the two CO(3-2) velocity components, 

^ , again benefiting from the magnification due to the foreground lens, modeling of which indicates that the offset 

^2 ' corresponds to just ~ 3 kpc in projection at z = 2.516. The spatial and velocity properties of these two components 

' are closely related to features detected in previously published Hq spectroscopy. We propose that this source is a 

^ . compact merger of two typical Lyman-break galaxies with a maximal separation between the two nuclei of about 

' 3 kpc, although a dusty disk explanation is not excluded. This system is much less luminous and massive than 

rS I other high-redshift submillimetre galaxies studied to date, but it bears a close similarity to similarly luminous, 

. dusty starburst resulting from lower-mass mergers in the local Universe. 
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1. Introduction 

The intensity of the far-infrared background indicates that 
far- infrared luminous, dusty sources make a significant 
contribution to the overall energy output from galaxies 
over the lifetime of the Universe (Puget et al. 1996, Fixsen 
et al. 1998). This emission may either represent AGN ac- 
tivity, in which case these dusty sources must coincide 
with an important phase in the growth of super-massive 
black holes, or dust-obscured star formation, in which case 



* Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by INSU/CNRS 
(France), MFC (Germany) and IGN (Spain). 



they are responsible for much of the massive star forma- 
tion and metal enrichment occurring in galaxies at high 
redshifts. To understand the nature and properties of this 
activity we need to resolve the background into individ- 
ual sources, suitable for more detailed study. For sources 
contributing to the peak of the far-infrared background, 
100-200^m, this work relies on space-based observatories, 
with modest apertures and correspondingly bright con- 
fusion limits, resulting in the resolution and identification 
of only a small (and hence perhaps unrepresentative) frac- 
tion of the galaxies contributing to the background, < 10% 
(Dole et al. 2001). 
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Redshift 

Obs.Freq. (GHz) 2.505 2.510 2.515 2.520 




LSR velocity (km/s) 



Fig. 1. (Left) The IRAM/PdBI i2C0(3-2) map of SMM J16359+6612 (=SMM1) superposed on the HSTF8UW image 
of A 2218. We clearly identify the images SMMl-A, B and C, all three well centered on their optical counterparts. 
The contours start at ±2ct and are plotted in steps of ct = 0.24 mJy. The large dashed circle give the half-power 
size of the primary beam of the PdBI at 98 GHz at the three pointing positions. The synthesized beam is displayed 
in the lower-left corner and has a FWHM of 5.9" x 5.4" and a position angle of 98°. We stress that the effective 
beam size in the source plane is roughly 20 times smaller along the shear direction (PA of ^ 40°) in the vicinity of 
SMMl-B. The inset at top right shows the spectral energy distribution of SMMl-B from the 3 mm, 1.3 mm, 850^m 
and 450/im measurements (the 450- and 850-^m fluxes density is measured from the SCUBA discovery maps in Kneib 
et al. 2004a). (Right) This panel displays (from top to bottom) the velocity proflle of the three different submm images 
SMMl-A, B and C, and of the sum of these three spectra. A double-peak proflle is clearly observed for SMMl-B and 
-C, as well as in the sum of all three components (bottom panel). Fitting two Gaussians to the the sum of the three 
velocity profiles indicates that the two components are separated by ~280kms~^. There is weak evidence that the 
higher-velocity component has a higher fiux density, but a narrower width. For comparison, we mark on the redshifts 
of the two features identified by Kneib et al. (2004a) from their analysis of the HST imaging and Ha spectroscopy of 
this system: a and /? (see Figure 2). 



Atmospheric windows at longer wavelengths allow 
larger-aperture ground-based telescopes to resolve a 
greater fraction of the cosmic background at mil- 
limetre (mm) and submillimetre (submm) wavelengths. 
Continuum surveys in the submm/mm wavebands using 
both SCUBA on the JCMT and MAMBO on the IRAM 
30-m have identified several hundred of submm/mm galax- 
ies (hereafter SMGs) to date. However, again the rela- 
tively coarse resolution of these telescopes, ^ 10", results 
in source confusion limiting these observations to identi- 
fying individual sources brighter than > 2 mJy at 850 /im 
(e.g. Blain et al. 2002). Thus most of the best-studied sam- 
ples of SMGs have 850-^m fluxes in the range 4-10 mJy 



(e.g. Ivison et al. 2002; Webb et al. 2003; Bertoldi et al. 
2000) , which have also been the focus of recent work to de- 
termine the redshift distribution of the bright submm pop- 
ulation, yielding a median redshift of z — 2.3 (Chapman 
et al. 2003, 2004). Unfortunately, once again, this popu- 
lation may not be completely representative as they only 
contribute around 30% of the far-infrared background at 
submm wavelengths (Small et al. 2002). 

To access the submm sources which are responsible 
for the bulk of the cosmic submillimetre background (and 
hence contain most of the star-formation or AGN ac- 
tivity) we need to probe down to flux densities below 
ImJy at 850 ^im (Small et al. 2002, Cowie et al. 2002). 
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These sources are fainter than the confusion Hmit of blank- 
field submm surveys making them impossible to iden- 
tify or study individually in detail from such samples. 
Nevertheless, they are amenable to study by exploiting 
the natural gravitational magnification of massive clus- 
ters of galaxies (Kneib et al. 1996; Blain 1997). The mag- 
nification provided by the cluster lens both boosts the 
brightness of the source and provides a finer effective beam 
size - reducing the effects of confusion - and bringing the 
sub-mJy submm population within reach of current tele- 
scopes. This technique was used for the first extragalac- 
tic submm survey (Small et al. 1997; Blain ct al. 1999; 
Chapman et al. 2000; Cowie et al. 2002) and has pro- 
vided a powerful route to understand the nature of this 
population of distant, faint galaxies (Blain et al. 2002) as 
well as other high-redshift galaxy populations (e.g. Baker 
et al. 2004). These studies show that SMGs with 850-/im 
fluxes within a factor of two of 1 mJy have a surface den- 
sity of ~ 3-4 arcmin"^ and produce roughly half of the 
submm background (Small et al. 2002). Combining these 
SMGs with the population of brighter sources accessible in 
blank-fleld surveys, we can trace the properties of SMGs 
that contribute the majority of the background radiation 
at submm wavelengths. 

The first redshifts of SMGs were obtained for sources 
lensed by massive galaxy clusters (Ivison et al. 1998, 2000) 
and cluster lenses continue to aid in the spectroscopic 
identification of this population (e.g. Fraycr ct al. 2003; 
Kneib et al. 2004a; Borys et al. 2004). These redshifts 
enable us to search for the corresponding line emission 
from cold gas within these galaxies - providing valuable 
information on the dynamics of these systems, the role 
played in these by the baryonic component and the extent 
of the gas reservoir available to power future activity. The 
first such observations targeted the rotational transitions 
of i^CO gas with the OVRO (Frayer et al. 1998, 1999, 
Ivison et al. 2001) and IRAM arrays (Downes & Solomon 
2003; Neri et al. 2003; Gcnzcl et al. 2003). Similar work 
on brighter blank-field SMGs has built upon the recent 
SMG redshift survey of Chapman et al. (2003, 2004), in 
a major programme to detect CO line emission from a 
large sample of bright submm galaxies with the IRAM 
PdBI (Neri et al. 2003; Greve et al. 2005). However, all 
of the submm galaxies studied in CO so far (both clus- 
ter lens and blank-field targets) are intrinsically luminous 
systems, none are from the mJy-submm population which 
dominates the far-infrared background. 

Recently, Kneib et al. (2004a) and Borys et al. (2004) 
have discovered two strongly magnified and multiply- 
imaged SMGs at z = 2.516 and z = 2.911 lensed respec- 
tively by the massive clusters A 2218 and MS 0451-05. 
These sources have intrinsic submm fiuxes below 1 mJy, 
but with magnifications greater than 10 for at least two of 
the three detected images in each system. Their CO emis- 
sion should be within reach of interferometric observations 
with current instruments. 

In this paper, we present recent observations with the 
IRAM array of redshifted molecular CO emission in the 3- 



mm and 1.3mm bands of the triply- imaged sub-mJy SMG 
SMMJ16359+6612 at ^ = 2.516. Interferometric CO ob- 
servations of this system have recently been published 
by Sheth et al. (2004) using the OVRO array, although 
the significance of their detections is lower than that de- 
scribed here. We present our observations in §2, describe 
our analysis and results in §3 and give our conclusions in 
§4. Throughout the paper we will assume an = 0.27, 
A = 0.73 cosmology with i?o = 71kms"^Mpc"^ At 
a redshift oi z = 2.516 the angular scale is thus 8.18 
kpc/arcsec. All quoted coordinates are in J2000. 



2. Observations 

We used the six antennas of the IRAM PdBI to image the 
CO(3-2) and CO(7-6) lines of SMM J16359+6612 (named 
hereafter SMMl in the text) at a redshift of 2.516 (Kneib 
et al. 2004a). As the primary beam of the antenna is 
about 50" and 21" at the frequency of redshifted CO(3-2) 
and CO(7-6) (98.405 GHz and 229.544 GHz respectively), 
three different positions around SMM J16359+6612 were 
targeted to completely map the three different images. We 
have approximately 16 hrs on-source at 3mm and 5 hrs at 
1.3mm (after fiagging low-quality visibilities) for the cen- 
tral pointing at a=16:35:54.5, (5=66:12:31, ~9hrs at 3mm 
on the southern pointing (q;=16:35:51.7, (5=66:12:10) and 
~ 2 hrs on the northern one (q=16:35:55, (5=66:12:50). 
No useful integration time was obtained at 1.3mm at the 
latter two positions. All these observations were carried 
out in the D configuration between November 2003 and 
June 2004. To improve the 1.3mm continuum detection, 
we also used data from an earlier ~5hrs track which we 
observed on December 10, 2001 at 237.191GHz with the 
phase tracking center at a=16:35:54.1, (5=66:12:23. We 
used 1637+574 as phase and amplitude calibrator, and 
0420-014, 0923+392, 3C84 and MWC349 as bandpass 
calibrators and to define the fiux density scale. We esti- 
mate that the fiux densities at 3mm are accurate to better 
than 10%, and to better than 20% at 1.3mm. 

Our 3-mm spectra provide velocity coverage extending 
from —850 to +850kms~^ around a zero-LSR redshift of 
z = 2.514 for each of the three sources. 3-mm line visibili- 
ties were averaged in frequency and combined to produce 
natural-weighted dirty maps. These were cleaned using 
the Clark algorithm, restored with a Gaussian beam of 
5.9" X 5.4" with a position angle of 98°, and corrected for 
primary beam attenuation. Natural-weighted dirty maps 
were also made of the continuum emission at 3mm and 
1.3mm by uv-averaging over the line-free channels. We 
estimate the absolute precision of the astrometry of our 
maps as 0.2-0.3", sufficient to precisely locate the CO 
emission relative to sources in the H.ST imaging. The as- 
trometry of the H.ST imaging is tied through wide-field 
ground-based images to the USNO reference frame, yield- 
ing absolute astrometry at the 0.2" level (Kneib et al. 
2004a). 
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Table 1. Observed properties of the three images of SMM J16359+6612, uncorrected for lens magnification. We hst 
the centroids of the CO(3-2) emission and the fiux-weighted redshifts derived for each image. /co(3-2) is the velocity- 
integrated fine intensity from the Gaussian fit (see text). Ib/Ir is the corresponding ratio of the blue-shifted to 
red-shifted CO(3-2) line emission. /co(7-6) /-^co(3-2) is the intensity ratio of CO(7-6) versus CO(3-2) line emission in 
the 0-370 km s~^ window (with respect to z = 2.514, Figure 1). S'l.smm is the fiux density of the 1.3mm continuum 
detected towards SMMl-B, the brightest component. The last column gives the magnification factor at the position 
of the three images derived from the detailed mass model of the foreground cluster by Kneib et al. (2004a, b). 

ID a,S Z -^CO{3-2) Ib/Ir -''cO(7-6)/^CO(3-2) Si. 3mm Mag 

(J2000) (Jykms"^) (mjy) 

A 16:35:55.18 66:12:37.2 (±0.2") 2.51740 1.67 (±0.13) 0.74 (±0.02) 2.8 (±1.2) ... 14 (±2) 
B 16:35:54.10 66:12:23.8 (±0.2") 2.51741 2.50 (±0.12) 0.74 (±0.02) 1.4 (±0.4) 3.0 (±0.7) 22 (±2) 
C 16:35:50.96 66:12:05.5 (±0.3") 2.51744 1.58 (±0.17) 0.88 (±0.02) 9 (±2) 




Fig. 2. 10"xlO" regions from the F814W-band HST WFPC2 image of the triply-imaged submm source 
SMM J16359±6612, showing images A, B and C (from left to right). We overlay on each the position of the CO(3-2) 
fine emission centroid as measured from our data (small circle), the relative uncertainty on these positions is i^ 0.3". 
For component SMMl-B (middle panel) we plot the size of the beam-corrected CO(3-2) profile as an ellipse and we 
also indicate the centroid of the blue-shifted (<220kms~^) and red-shifted (>280kms^^) velocity components of the 
line (open squares) and the location of the 1.3mm continuum emission (large open cross). Note that every image of 
SMMl comprises a NIR source (7, marked as a cross) which is bracketed by two bluer features visible in these F814W 
images (a and /3). For SMMl-A and SMMl-B we have indicated by white lines the extent of the Ha fine observed in 
the NIRSPEC data of Kneib et al. (2004a). The resolution of the HST images are '--0.15"; North is up and East is to 
the left. 



3. Analysis and Results 

3.1. The CO (3-2) line 

We detect the CO(3-2) line at high significance in all 
three images and show the velocity-integrated line map 
and spectra in Figure 1, with tabulated information in 
Table 1. The fiux-weighted mean redshifts of the fines (in- 
tegrated across the FWZI from —10 to 540km s~^) in the 
three images are identical - zco = 2.5174±0.0002. The ra- 
tio of CO(3-2) line intensities for the three images roughly 
corresponds to that seen for their 850-/xm and K-hand 
fiuxes from Kneib et al. (2004a): Ia/b = 0.67 ± 0.18 and 
/^./B = 0.63 ±0.19, compared to 0.64, 0.53 and 0.53, 0.40 
from the 850-/im and K-hand respectively. These fiux ra- 
tios are all consistent with the predicted magnification ra- 
tios from the lens model of the cluster: Ma/b = 0.64±0.14 
and Mc/B = 0.41 ± 0.12. The CO(3-2) line widths mea- 
sured from the spectrum of each image are also all large. 



corresponding to FWZI of 500 km s"^. These observations 
allow us for the first time to categorically state that all 
three sources, SMMl-A, B and C, represent images of a 
single background object. 

Compared to the OVRO observations of 
SMMJ16359±6612 in Sheth et al. (2004), we mea- 
sure a slightly lower total integrated CO fiux: 
5.75 ± 0.25Jykms-i (versus 6.3 ± 0.2Jykms-i), 
and a slightly higher redshift: z = 2.5174 ± 0.0002 
(compared to z = 2.5168 ± 0.0003). In both cases the 
differences are relatively minor and the higher signal to 
noise of our detections would suggest our results are the 
more reliable. 

We can also exploit the better signal-to-noise of our 
observations to look in more detail at the structure in 
the CO (3-2) fines. We find that the CO (3-2) fine profiles 
display a characteristic double-peaked profile in SMMl- 
B and SMMl-C, while the structure of the spectrum of 
SMMl-A is less obviously bimodal. We combine the three 
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Table 2. The properties of the velocity components of the CO (3-2) and Ha lines of SMMl-B. The coordinate offsets 
are computed with respect to the position of the optical feature a within SMMl-B (q:=16:35:54.19, 5=66:12:24.92, 
Kneib et al. 2004a) - note there the relative astrometry of the optical and CO maps is uncertain at the ^ 0.4" level. 
CO(3-2) and Ha velocities are relative to the Ha redshift oi a oi z = 2.5165 (corresponding to +210kms~^ relative 
to .2 = 2.5140). 



SMMl-B 


Aa 


A6 


2 


V 


Av 


Component 


(") 


(") 




km s~^ 


kms^^ 


CO(3-2) blue+red 


-0.56 (±0.16) 


-1.14 (±0.16) 


2.5172 (±0.0002) 


70 (±20) 


440 (±80) 


CO(3-2) blue 


0.42 (±0.20) 


-0.37 (±0.20) 


2.5156 (±0.0002) 


-70 (±20) 


220 (±50) 


a (Ha) 


0.0 


0.0 


2.5165 (±0.0015) 





280 (±60) 


CO(3-2) red 


0.45 (±0.17) 


-1.41 (±0.17) 


2.5189 (±0.0002) 


210 (±20) 


220 (±50) 


13 {Ha) 


0.78 (±0.10) 


-2.51 (±0.10) 


2.5190 (±0.0015) 


220 (±60) 


280 (±60) 



spectra to increase the signal-to-noise yet further (as- 
suming that differential magnification effects across the 
source are small, which should be correct for a compact 
object) as shown in Figure 1. We fit the combined ve- 
locity profile using a two Gaussian fit to determine the 
redshift, width and relative intensity of the two veloc- 
ity components. We find that the velocity profile requires 
both Gaussians for an adequate fit, each with widths of 
220±20kms-i and separated by 280±20kms"^ The fiux 
ratio between the blue-shifted and red-shifted components 
is Ib/Ir = 0.74 ± 0.02 (Table 1). The redshifts of the 
two velocitjf components arc Zbiuc = 2.5156 ± 0.0002 and 
Zred = 2.5189 ±0.0002. We note that Kneib et al. (2004a) 
reported z = 2.5165 ± 0.0015 for the Ha redshift of the 
optical structure a, and z = 2.5190 ± 0.0015 for feature ,,9 
in SMMl-B (Figure 2) and return to this point below. 

To investigate these two velocity components in more 
detail and search for changes in size or position with 
velocity, the visibilities were fitted first with elliptical 
Gaussians. Although we could not resolve SMMl-A or -C, 
SMMl-B (the most magnified image) is best fit by an el- 
liptical Gaussian with a beam-corrected size of 4.1 ±0.6" x 
1.8 ± 1.0"and a position angle of 30 ± 13° - this is plot- 
ted as an eUipse in the central panel of Figure 2. We note 
that this agrees well with the size and orientation of the 
optical counterpart of this system in Figure 2. Correcting 
for the lensing magnification, using the well-constrained 
mass model of A 2218 (Kneib et al. 1996, 2004b; EUis et 
al. 2001), this ellipse maps to a linear size of ~ 3 x 1.5 kpc 
in the source plane (there is some complex, but weak, de- 
pendence according to the adopted position angle). Note 
that the tentative resolution of SMMl-B claimed by Sheth 
et al. (2004) has a nearly orthogonal orientation to that 
derived here (PA of —68°), transverse to the local shear 
direction and hence the lensing-corrected size they esti- 
mate is significantly larger than our measurement: 12 kpc 
versus 3 kpc. 

To investigate whether the red-shifted and blue-shifted 
velocity components originate from two different spatial 
positions, for the blue-shifted emission we integrated the 
velocity map for all three images from ~ 20-230 km s~^ 
(note the reference redshift is 2; = 2.5140) and for the 



red-shifted emission from ~ 300-510 km s ^. Again for 
SMMl-A and SMMl-C, no strong distinction could be 
made in the position of the blue-shifted and red-shifted 
velocity channels. However, for the SMMl-B component 
the two velocity components show a well-detected spa- 
tial separation (see Table 2 and Figure 2) between the 
blue-shifted and red-shifted gas. The two components are 
separated along the major axes of the CO (3-2) emission 
and the optical counterpart (Figure 2). Based on just the 
spatial information, and allowing for the ~ 0.4" uncer- 
tainty in the relative astrometry of the CO and optical 
maps, we suggest that the blue-shifted component of the 
SMMl-B CO(3-2) fine is Hkely associated with the optical 
feature a and the red-shifted component with the opti- 
cal feature (3. If this is indeed the case, then the blue- 
and red-shifted components of the CO (3-2) emission in 
SMMl-B are both spatially and kinematically coincident 
with the optical structures a and /? cataloged by Kneib et 
al. (2004a). Alternatively, the CO(3-2) line could be asso- 
ciated to a smaller rotating structure closer and around 
component 7, however with no strong optical or Ha coun- 
terpart, which makes this possibility less likely but not im- 
possible. Higher resolution millimeter observation of this 
system should confirm the true nature of this system. 

3.2. The CO(7-6) line 

The velocity coverage of our parallel observations of the 
CO(7-6) SMMJ16359+6612 in the 1.3mm band extends 
from —380 km s^^ to +380kms~^ (corresponding to the 
redshift interval z = 2.510-2.518). Thus, only part (about 
70%) of the expected broad CO(7-6) line is covered by 
the frequency setup. The noise level in the 1.3mm band is 
also relatively high and does not allow a high confidence 
level measurement. Within this limited 1.3mm window we 
measure a CO(7-6) fiux of /co(7-6) = 3.3 ± 1.4 Jy.km/s 
for SMMl-A and /co(7-6) = 2.5 ±0.7 Jy.km/s for SMMl- 
B. In order to quote representative measurements of the 
CO(7-6)/CO(3-2) line ratio, we have computed the inte- 
grated flux in the redshift window from z = 2.5140 to 
z = 2.5185 for the two lines CO(3-2) and CO(7-6) at the 
position of SMMl-A and SMMl-B. The derived flux ra- 
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tios are: /co(7-6)/-fco(3-2) = 2.8 ± 1.2 for SMMl-A, and 
-fco(7-6)Aco(3-2) = 1.4 ± 0.6 for SMMl-B (the larger 
error for SMMl-A is mainly due to the larger primary 
beam correction at 1.3mm for this source). These mea- 
surements should be reduced by ~ 50% to account for the 
contribution from continuum emission 3 mJy) within 
the velocity window used to integrate the CO (7-6) line. 

The mean CO(7-6) to CO(3-2) flux ratio for SMMl- 
A/'B is comparable to the value /co(7-6)/^CO(3-2) = 
1.14 ±0.20 found for the Icnscd SMG, SMM J14011+0252 
at z = 2.565 (Downes & Solomon 2003). As in 
SMMJ14011+0252, the strength of CO(7-6) relative to 
CO (3-2) indicates both a relatively high temperature, 
T ~ 50 K, and high H2 density, ~ 10"^ cm^'\ for the gas 
reservoir in this galaxy. 

3.3. Continuum emission 

We detected weak continuum emission towards SMMl- 
B at a =16:35:54.1, S =66:12:24.36. The flux density is 
3.0 ± 0.7 mJy at 1.3mm (237 GHz), and within our as- 
trometric accuracy (0.2"-0.3") this position (Figure 2) 
is coincident with the integrated CO(3-2) emission from 
SMMl-B (the synthetised beam has a size of 1.9" x 1.5" at 
a PA of —37°). We also detect a much weaker continuum 
source at this precise position in our 3 mm (98 GHz) data, 
with a flux density of 0.29 ± 0.13 mJy. We flt the spec- 
tral energy distribution (SED) of the dust emission using 
these two flux points, in combination with the 450- and 
850-/im values from Kneib et al. (2004a), see inset panel in 
Figure 1. We determine a spectral index of ~ v^'^ from the 
slope of the SED between 850 /im and 3 mm - character- 
istic of optically thin dust emission. The dust emissivity 
index is thus (3 ~ 1.2, very similar to that derived for 
bright IRAS galaxies locafly {(3 ~ 1.3 ± 0.2, Dunne et al. 
2000). 

We can also estimate the temperature from the fit 
to the four long-wavelength photometric observations of 
the SMMl-B, this is determined rather accurately, as 
Td = 51 ± 3K. Including the Metcalfe et al.(2003) ISO 
hmit we have = 48 ± 3K for SMMl-B data assuming 
(3 = 1.5. We can use this value, along with the dust emis- 
sivity index and the magnification-corrected 850 /xm flux 
(which is our most precise flux measurement) to estimate 
the dust mass, assuming a dust mass opacity coefiicient, 
Kd = 0.077 m^Kg-i at 850 /xm (Dunne et al. 2000). This 
calculation yields a mass of (1.9 ± 0.3) x lO'' Mq, subject 
to the considerable uncertainty in the opacity coefficient. 

The corresponding total bolometric luminosity is 1.6 ± 
0.4 X 10^^ Lq. Hence, although intrinsically faint at mil- 
limetre and submillimetre wavelengths, owing to its high 
dust temperature, SMM J16359+6612 is still a relatively 
luminous galaxy. It will be interesting to include Spitzer 
Space Telescope data to confirm and better constrain the 
dust temperature but unfortunately there is currently no 
deep radio image of A 2218. 



4. Conclusions 

We have compared the results of our PdBI observations 
of SMM J16359+6612 with those derived from shaflower 
data taken with the OVRO array by Slicth ct al. (2004). 
We find broad agreement between these datasets, some 
discrepancies do exist but can be attributed to the lower 
signal to noise of the OVRO data. The high-fidelity CO 
maps provided by the IRAM PdBI (in particular due to 
its better coverage of all 3 sources) add significantly more 
detail to address the nature and properties of this unique 
source. 

A detailed analysis of the CO(3-2) velocity profile 
shows a double-peak structure with a significant spatial 
offset between the blue- and red-shifted gas components. 
Comparing the positions and redshifts of these two gas 
components with the complex optical /'near-infrared mor- 
phology discussed by Kneib et al. (2004a), we tentatively 
identify the kinematic components seen in CO(3-2) with 
two features visible in the HST imaging of this system: 
a and /3 (Kneib et al. 2004a). We note that the C0(3- 
2)/Ha luminosity ratio is much larger for component f3, 
either due to stronger dust obscuration in or because 
of slit-losses in the Ha observations resulting from a mis- 
alignment of the slit across the source. Near-infrared inte- 
gral field spectroscopy of this SMG to measure the spatial 
structure of the Ha fine wfll be essential to better de- 
termine the relative distribution of cold and ionised gas 
across this system. 

The spatial separation of the two CO(3-2) components 
is ~ 1.5 kpc (with a maximal extent of ~ 3 kpc) on the sky 
(once corrected for the lensing magnification). This scale 
suggests we are either seeing a molecular gas disk or ring 
within a single galaxy, or two strongly interacting gas-rich 
components of a merger. The latter possibility is favored, 
as the three optical/NIR component a,/3 and 7 are not 
co-hnear, and thus are not representative of a stable disk. 
We can use the spatial and velocity information (a rotation 
velocity of 150kms~^ at a scale of 3 kpc) to estimate the 
dynamical masses in these two scenarios. We derive a mass 
of Mdyn = (0.8 ± 0.2) X IQIOM© at a radius of 1.5 kpc if 
the gas is distributed in a ring/disk, and M^yn = (1.5 ± 
0.3) X 10^° M0 within the same aperture if the system is a 
merger (we assume a 90° inchnation here, otherwise both 
of these estimates need to be scale with (sini)~^). 

SMM J16359+6612 is the third SMG in which 
spatially-resolved CO emission has been detected, the 
other two being SMM J02399-0135 at z = 2.81 (Genzel 
et al. 2003) and SMM .J14011+0252 at z = 2.56 (Downes 
& Solomon 2003; Ivison ct al. 2001). SMM J02399-0135 
shows a double-peaked CO (3-2) fine with a spatial offset 
of 10 kpc and 400kms~^ between the components, while 
SMM J14011+0252 has only a single velocity component 
(FWHM of 200kms-i) with a size of ~ 3 kpc (adopting 
the magnification factor for the lensing configuration for 
this source from Swinbank et al. 2004). The gas distribu- 
tion in SMM J14011+0252 is thus very similar in scale (as 
well as having similar dynamical and far-infrared lumi- 
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nosities) to SMM J16359+6612, while SMM J02399-0135 
is both significantly more luminous, larger and more mas- 
sive. However, the dense gas detected in all three sys- 
tems has a much wider distribution than the compara- 
ble component in local Ultraluminous Infrared Galaxies 
(ULIRGs), - 0.3 kpc (Downes & Solomon 1998). 

With our assumed cosmology, the total CO (3-2) line 
flux converts to an observed CO(3-2) luminosity of SMMl- 
B of (8.2 ± 0.4) X 10-^" K kms~^ pc^, correcting from the 
lensing magnification the intrinsic CO(3-2) luminosity is 
then (3.7 ±0.2) x 10^ K km s'^ pc^. The CO(3-2) lumi- 
nosity we measure is roughly half of the value found for 
Arp 220 (Solomon et al. 1997), although Arp 220 is compa- 
rably luminous in the far-infrared to SMM J16359+6612. 

Using the relation from Solomon et al. (1997) and 
Downes & Solomon (1998) which indicate a CO lumi- 
nosity to gas mass conversion factor of ^ 0.8-1.6 M© 
(Kkms""'^ pc^)~^ we derive a molecular gas mass of (4.5± 
1.0) X 10^ Mq, giving a gas mass surface density of > 
1O^M0 pc~^. Comparing this to our dynamical estimates 
suggests a gas mass fraction of ^ 60% for a disk /ring- like 
structure or ~ 30% for a merger. Combined with our ear- 
lier estimate of the dust mass, this indicates a dust-to-gas 
ratio of 0.004 ±0.001, similar to that seen in local galaxies 
(Dunne et al. 2000). 

The star-formation rate estimated from the 
far-infrared luminosity by Kneib et al. (2004a) is 
~5OOM0yr~^ for stars >O.1M0 (note that in contrast 
the Ha-derived SFR is only ~ 10 Mq yr^^, before correct- 
ing for extinction and slit losses). Taking the far- infrared 
estimate of the SFR this indicates that the current 
activity cannot continue at its present level for more 
than another lOMyrs. Even allowing for the fact that the 
far-infrared luminosity is dominated by the emission from 
the most massive stars - thus allowing us to use an IMF 
dominated by massive stars or even to adopt a top-heavy 
IMF (Baugh et al. 2005) - would only extend this lifetime 
to ~ 50 Myrs. These estimates are comparable to the 
crossing-time of the system, indicating that the starburst 
is likely to be an instantaneous non-equilibrium event - 
arguing against a stable disk/ring configuration for the 
gas within this system - and suggesting instead that 
we are seeing the results of a tidal-induced starburst in 
a pair of merging galaxies (a and /3). The rest-frame 
UV/optical morphology of this system suggests that a 
highly-obscured starburst is being triggered in the overlap 
region between the two galaxies (component 7 in Kneib 
et al. 2004a) as is commonly seen in low-redshift mergers. 

We can also compare the ratio of the far-infrared lu- 
minosity produced by SMM J16359+6612 to its H2 gas 
mass: Lfir/-Mh2 ~ 320. This is roughly ten time higher 
than seen in Arp 220 (or comparably luminous galaxies lo- 
cally, Dunne et al. 2000; Downes et al. 1997). This high 
ratio indicates that the activity in SMM J16359+6612 is 
far more efficiently or vigorously using the available gas 
supply than is typical in starburst systems at 2 ~ 0. This 
may either result from special physical conditions associ- 
ated with intense starbursts at high redshifts. For example 



the activity could derive from wide spread instabilities in 
large gas disks triggered by dynamical interactions, which 
are not possible in comparably massive and gas-rich, but 
bulge-stabilised, systems at low redshift. Alternatively, the 
high Z/fir/-^H2 could reflect far-infrared emission from 
a highly-obscured AGN - although there arc no spectral 
signatures of an AGN. Turning this argument around - 
the free availability of gas to power such nuclear activity, 
but the apparent lack of any signatures of an AGN in the 
spectroscopy of this source from Kneib et al. (2004a), may 
indicate that the central black hole(s) in this system (and 
the associated bulge) has not yet grown large enough to 
produce a luminous AGN. The current interaction, and 
any subsequent mergers with the other close compan- 
ions in the compact and dynamically-cold group around 
SMM J16359+6612 (the Lyman-break galaxies #384 and 
#273, Kneib et al. 2004a), will provide the opportunity for 
substantial bulge growth and associated accretion-driven 
activity from the central black hole of this system. 

To judge how the properties of this relatively low- 
luminosity SMG compare to the other well-studied popu- 
lation of moderately active galaxies aX z > 2, the Lyman- 
break population (Steidel et al. 2004), we can compare the 
parameters we measure for SMM J16359+6612 to those 
measured for a similarly highly magnified Lyman-break 
galaxy: cB58 a.t z = 2.73 (Baker et al. 2001, 2004). cB58 
has a SFR of just 2OM0yr~^, an inferred gas mass of 

7 X 1O9M0, a dynamical mass of ~ l-2xlOi"M0, a 
gas fraction of ~ 0.3-0.7 and an average gas surface mass 
density of S ~ 5OOM0pc~^ (all within a 2 kpc radius. 
Baker et al. 2004). These estimates are very similar to 
those we derive for SMM J16359+6612 with the exception 
that SMM J16359+6612 has an order of magnitude higher 
instantaneous star-formation rate and a marginally higher 
gas surface density. Taken with our previous suggestion 
that SMM J16359+6612 represent a merger, the compari- 
son to cB58 suggests that it may be a merger between two 
typical Lyman-break galaxies. The merger has increased 
the gas density within the system and triggered a burst 
of activity which will quickly use up the bulk of the gas 
reservoir in a single event (which would have otherwise 
fueled a more sedate level of activity for lOO's Myrs). 

In summary, several of the properties of 
SMM J16359+6612 appear to resemble those of sim- 
ilar luminosity merging/interacting systems at low 
redshift, e.g. the dust spectral index, dust-to-gas ratio, 
presence of an obscured starburst in the overlap zone 
between the merging components, etc. However, the ratio 
of far-infrared luminosity to gas mass and the physical 
extent of the apparently high-densitj^ gas reservoirs in 
the progenitors are both an order of magnitude larger 
than for local systems. As SMM J16359+6612 is an order 
of magnitude fainter at submm wavelengths than other 
distant SMGs studied in detail to date, and also less 
massive, it is not yet clear how many of these similarities 
and differences to low-redshift ULIRGs will carry over 
to the bulk of the SMG population. However, we stress 
that the properties of SMM J16359+6612 may be more 
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representative of the high-redshift population which pro- 
duces a large fraction of the submm cosmic extragalactic 
background in the submm waveband, and hence are more 
relevant for understanding the evolutionary behaviour 
of the majority of the oldest stellar components in local 
galaxies: bulges. The fact that the properties of this 
system are more likely to resemble those expected for a 
merger of two typical Lyman-break galaxies suggests a 
Hnk between these two populations (at least for sub-mJy 
SMGs). We also note that SMM J16359+6612 has yet to 
be observed at the highest spatial resolution achievable 
with the IRAM array - these will provide a further factor 
of ten improvement in the resolution - probing the gas 
distribution within this rare example of the sub-mJy 
submm galaxy population at sub-kpc scales, which will 
fully characterize the nature of this system. 
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